Charge order versus Zener polarons: ferroelectricity in manganites. 
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We show that in manganites close to half-doping novel non-bipartite magnetic phases appear 
due to the interplay between double exchange, superexchange and orbital ordering. In considerable 
part of the phase diagram the groundstate has a magnetic order that is intermediate between the 
canonical magnetic CE-phase and a state that we identify as the recently observed Zener polaron 
state. The intermediate phase shows a type of charge ordering that breaks inversion symmetry and 
is therefore predicted to be ferroelectric. 
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Introduction. The wide variety of phenomena that are 
observed in transition metal (TM) oxides is caused by 
the interplay of various degrees of freedom. An intrinsic 
property of doped TM oxides is the formation of super- 
structures that are related to the charge, spin, lattice and 
orbital degrees of freedom. Different types of charge or- 
dering (CO), orbital ordering and magnetic ordering may 
exist depending on e.g. doping concentration, tempera- 
ture and pressure. Charge order is usually considered 
to be an ordering of TM ions with different valence, for 
instance the alternation of Mn 3+ / 4+ in Lao.sCao.sMnOa 
or Fe 2+ / 3+ in Fe304. Such a CO state we refer to as 
a site-centered (or TM-centered) CO, or equivalently, a 
site-centered charge density wave (SCDW). There is how- 
ever another possibility for charge ordering — the forma- 
tion of a bond-centered charge density wave (BCDW). 
The BCDW is very similar to the state that is observed 
in Peierls distorted quasi one-dimensional materials. In 
a TM oxide the BCDW is actually an oxygen-centered 
charge ordered state as oxygen ions are located on the 
TM-TM bonds in the typical (e.g. perovskite) crystal 
structures of these systems. The possibility to have such 
an oxygen-centered BCDW was studied theoretically for 
the compounds RNi0 3 (with R =Pr, Nd, etc.) 0|. 

For about the last 50 years it was presumed that in 
manganites of the type Ri_ x Ca x Mn03 (with R = La, 
Pr, etc.) only the Mn-centered SCDW occurs. At e.g. 
x = 1/2 this SCDW state has the well-known orbital 
ordered CE-type magnetic structure, see Fig. la, an or- 
dering for which there is abundant experimental and 
theoretical support 0, 0- However, Daoud-Aladine et 
al. challenge this canonical picture in their recent study 
0. These authors conclude that in single crystalline 
Pro.6Cao.4Mn03 another type of superstructure is more 
consistent with the experimental data. This superstruc- 
ture is a BCDW, which they refer to as the Zener polaron 
state (Fig. lb). This state is made out of dimcrs in which 
all Mn ions have the same valence (all Mn 3+I ), quite 
different from the canonical SCDW that holds Mn ions 
with unequal valence. Recently the BCDW state was 
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FIG. 1: Three types of magnetic/orbital superstructures: a) 
CE phase; b) Orthogonal (_L) phase; c) 120° - phase. Direc- 
tions of spins are shown by arrows. 



also found theoretically in Hartee-Fock calculations of 
the electronic structure 0,0. The presence of a BCDW 
amounts to a drastic change in our understanding of 
charge ordering — not only for manganites. It may imply 
that in all cases of charge order in oxides one should have 
in mind, besides the usual site-centered CO, alternative 
states with bond-centered superstructures. 

In this paper we consider the regions of stability of dif- 
ferent superstructures in doped manganites close to half- 
doping, and show that indeed the bond-centered charge 
density wave (Zener polaron) state may be favorable in 
part of the phase diagram. Our results reconcile the con- 
flicting conclusions of Rcfs. and 0]: the conventional 
CE state is still better at exactly half-doping x = 1/2, but 
a BCDW is stable at around x ~ 0.4, as observed exper- 
imentally. The most interesting result of our calculation 
is the following: in a large part of the phase diagram an 
intermediate state actually turns out to be most favor- 
able. This state is a superposition of a BCDW (Zener 
polaron) and a SCDW (of CE-type). As one can easily 
see, a superposition of Figs. la and lb results in a state 
in which the two Mn ions in one dimer are inequivalent. 
This produces a dipole moment on each pair, and the 
resulting state is therefore ferroelectric. The appearance 
of ferroelectricity in the case of coexistence of a bond- 
centered CDW (Peierls state) and the site-centered CDW 
was previously observed in some quasi one-dimensional 



compounds [jj. But, to the best of our knowledge, such 
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a state was never discussed in the context of the three- 
dimensional TM oxides. In this respect it is noteworthy 
that first measurements of the dielectric constant in these 
manganites show a strong anomaly at the CO phase tran- 
sition, which might be related to ferroelectricity [8j. 

Qualitative considerations. It is well established that 
charge and orbital structures are closely related to the 
magnetic ordering [|J. In all the papers where either the 
usual SCDW or the BCDW (Zener polarons) are 
considered, the magnetic structure is taken to be of the 
CE type (ferromagnetic zig-zag chains stacked antiferro- 
magnetically, Fig. la). However in the BCDW picture the 
situation can change drastically due to orbital ordering 
existing in the structure of Fig. lb (predominantly d x i- 
orbitals on horizontal dimers and d y 2 on vertical ones, see 
below) . Due to the double exchange mechanism the cou- 
pling between spins inside a Zener polaron is definitely 
strongly ferromagnetic. The high-spin polaron dimers — 
if we treat them as localized spins — form a triangular 
lattice, see Fig.lc, with a relatively strong antiferromag- 
netic coupling between neighboring parallel dimers. A 
different magnetic coupling can be expected when the 
dimers are are oriented perpendicular to each other. If 
these couplings would all be comparable in magnitude 
and antiferromagnetic, then the angle between neighbor- 
ing high spin polarons on the triangular lattice would 
be 120° and a tripartite Jaffet-Kittel magnetic structure 
would be formed, Fig.lc. However from the orbital struc- 
ture of Zener polaron described above one can immedi- 
ately conclude that the perpendicular coupling is defi- 
nitely weaker. We therefore expect that the magnetic 
structure tends to an "orthogonal" (_L) ordered state 
with four different spin sublattices as shown in Fig. lb. 
We should mention that the neutron scattering data for 
Pr .6Cao.4Mn03 can be equally well described both by 
the conventional magnetic structure of Fig. la and by 
the orthogonal structure of Fig. lb ^3 that we propose 
here. In the overdoped manganites Ri-jjCa^MnOa, with 
x = 2/3, 3/4 a _L magnetic structure was already ob- 
served [ii| . 

Calculated phase diagram. In our calculations we use 
a tight binding approach for the band structure of the 
manganites, where both double exchange and superex- 
change are incorporated. This approach proved to be 
quite efficient for overdoped (x > 0.5) [T3I and half- 
doped manganites. In the double-exchange frame- 
work with strong on-site Hund's rule coupling the mo- 
tion of e g -electrons is largely determined by the underly- 
ing magnetic structure. In this approach the stability of 
various magnetic and orbital structures is determined by 
the competition of the band energy of the e s -electrons, 
favoring ferromagnetism, and superexchange interaction 
J between localized (t2 g ) spins, favoring antiferromag- 
netism. The double degeneracy of conduction electrons 
significantly modifies the conventional picture and leads 
to the stabilization of more complicated magnetic struc- 
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FIG. 2: Phase diagram of manganites for doping concentra- 
tions close to 1/2. F: ferromagnetic, G: antiferromagnetic 
(AF) A: F planes coupled AFM, C: F chains coupled AF, 
CE: F zig-zag chains coupled AF, 120: Jaffet-Kittel state 
and I: intermediate orthogonal phase. The electron density is 
n e i = 1 — x. 

tures, besides the simple ferromagnetic (F) or antiferro- 
magnetic (G) ones 12] . We calculate the groundstate 
energy for different magnetic structures: G-, C- and A- 
and CE-type antiferromagnetic states, the ferromagnetic 
(F) phase, and, in addition, for two new phases (with 
orthogonal (_L) and 120° Jaffet-Kittel type magnetic or- 
dering of spin dimers) that we expect to be relevant close 
to x ~ 1/2, based on the arguments above. 

The effective hopping matrix elements tf? depend on 
the type of the orbitals a, (5 and relative orientation of 

2 2 2 2 2 2 

a pair i, j, so that ' x ~ v = 3/4i, i*^* = t and 

if. ij* = l/4i 0|- Besides this, there is the usual de- 
pendence of hopping on the angle 9 between neighboring 
spins: tf'^ = i™ ,/3 cos %/2. The antiferromagnetic su- 
perexchange energy per bond is J cos 6ij. With these 
expressions it is straightforward to calculate the band 
dispersion and total energy of different phases, cf. [l^ . 
In Fig. 2. we show the resulting phase diagram as a func- 
tion of doping x and antiferromagnetic coupling J. For 
x > 1/2 we find mainly the magnetic CE, C, A and F 
phase to be stable, in accordance with previous results 
[12|,[l5j . For a; < 0.5 we also find a large region where the 
BCDW magnetic 120° state is stable. Most interestingly, 
we also find a finite region of the _L phase. This phase 
is stable for doping concentration from around x ~ 0.4, 
which agrees with the fact that at this concentration the 
Zener polaron state was found 

From the calculation we know the corresponding wave 
functions so we can find the orbital occupation on a cer- 
tain site by projecting the wave function on that par- 
ticular site. In agreement with the standard picture we 
obtain that for the CE-structure the edge sites A, A' in 
Fig. la. are occupied by orbitals d x i and d y 2 , respectively. 
On the corner sites predominantly the d x 2_ y 2 orbital is 
occupied with about 10 % admixture of d z 2 density. 

Let us denote the general orbital state as \<j>) = 
coscj)/2\z 2 ) + smcf>/2\x 2 — y 2 ). For an isolated dimer that 
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points along the x direction (the equivalent D, D' sites 
in Fig.l.b) we then expect (j> = 2ir/3 for both sites in 
the dimer, which corresponds to occupied d x 2 orbitals on 
both sites. For an isolated dimer along the y direction 
(C, C sites in Fig.l.b.) we then expect cf) = —2n/3, 
corresponding to occupied d y 2 orbitals. From the cal- 
culation of the wave functions for the _L phase, where 
there is an interaction between neighboring dimers, we 
find for the x-dimer that <j> ~ 130°, which is mainly a 
ci^-state, as we had for a non-interaction dimer, but 
with some admixture of d x 2_ y 2. Correspondingly, for the 
y-dimers we find <f> ~ —130°, which amounts to a pre- 
dominantly d y 2 -state with some d x 2_ y 2 character. This 
admixture is important as it is driving the stability of 
the _L phase. The superexchange tends to align the spins 
of neighboring x-and y-dimers antiparallel. But in the _L 
phase these spins are only orthogonal, so that the hop- 
ping between the dimers is possible, although reduced 
by a factor l/\/2 from its maximum value. From this 
hopping between neighboring dimers, which leads to the 
admixture of d x i_ y 2 states into the groundstate wave- 
function, the system gains kinetic energy. Note that the 
hopping amplitude between d x 2 orbitals of neighboring x- 
dimers (or between d y 2 orbitals of neighboring y-dimers, 
respectively) is much smaller, so that the spins of neigh- 
boring dimers of the same kind tend to align themselves 
fully antiparallel. Thus we conclude that the stability 
of the _L state is caused by its orbital order that allows 
for an optimal compromise between electron derealiza- 
tion, causing ferromagnetic double exchange, and anti- 
ferromagnetic superexchange interactions. 

Charge ordering and ferroelectricity. In the CE-phase 
the edge and the corner sites are crystalographically in- 
equivalent. This results in different total electron oc- 
cupation numbers on edge and corner sites. Already 
when the on-site Coulomb interaction U between elec- 
trons in different orbitals but on the same site is included, 
the common picture of a checkerboard SCDW arises in 
this system, with a valence Mn 3 5+Aq of the corner sites, 
Mn 3.5-A 9 of the edge sitcs and Aq < o.l 13]. The CE- 

structure has a center of inversion symmetry and, conse- 
quently, there are no electrical dipoles moments present 
in the groundstate. This is also the case in the pure 
_L structure as both Mn ions in each dimer are equiv- 
alent (lq|. But there may exist an intermediate state 
that does have dipole moments: this intermediate state 
is a superposition of the CE and _L phase. This state 
has a magnetic structure in which the spins of x-dimers 
(spins within each dimer are ferromagnetic) are antipar- 
allel to one another, and the y-dimers, also antiparallel 
between themselves, are rotated with a certain angle 9 
with respect to the x-dimers. The angle 9 characterizes 
intermediate phase as it determines the amount of mix- 
ing between CE and _L that is present in the groundstate: 
9 = corresponds to the CE-structure, while 6 = tt/2 
corresponds to the _L one. 




FIG. 3: Charge disproportionation Aq and mixing angle <f> 
(see text) for the U = t in Hartree-Fock approximation in the 
intermediate ± and CE phases. 



We find that the region of stability of this intermediate 
phase is larger than that of a pure _L-phase, and extends 
from the concentrations x ~ 0.4 to x = 0.5; Fig. 2 shows 
for which values of x and J/t this mixed intermediate 
phase is stable. The angle # = 0forx = l/2, and it 
increases with decreasing x, reaching the value 9 = tt/2 
(_L phase) at x ~ 0.4, see Fig. 3. The intermediate phase 
is therefore in general a superposition of the CE and _L 
phase, except near the phase boundaries in Fig. 2. where 
the pure phases are stable. In the intermediate phase 
there exist simultaneously a dimer (Zener polaron) bond- 
centered CDW and a checkerboard site-centered CDW. 
This implies that in the intermediate state the valence of 
the two Mn's in a dimer is different, inversion symmetry 
is broken and each dimer attains a dipole moment. In 
this phase the dipole moments of the x-dimers point e.g. 
to the right in Fig.l and all dipole moments of y-dimers 
point e.g. up (depending on how precisely the inversion 
symmetry of the system is broken globally) causing each 
xy-plane to acquire a net ferroelectric moment. 

We calculate the charge disproportionation Aq due to 
an on-site Coulomb interaction U between the electrons 
in the Hartree-Fock approximation. In Fig. 3. Aq is 
shown as a function of doping. Charge ordering starts 
at x ~ 0.4 where the intermediate phase appears, and 
has its maximum at x = 1/2, which coincides with the 
transition to a pure CE-phase. 

It is experimentally observed that the CDW deforma- 
tions in neighboring planes are in phase with each other. 
This implies that there is a net polarization of the electric 
dipole moments when the system is in the intermediate 
phase, which leads to a ferroelectric ground state. Of 
course, in real materials ferroelectric domains may ap- 
pear. As also the conductivity of these materials is not 
so low, it may not be easy to observe the predicted ferro- 
electricity in CO manganites. But already first measure- 
ments of the dielectric constant e do show a strong peak 
in e(T) at the charge ordering temperature 0. A more 
detailed experimental study of this phenomenon would 
certainly be of interest in the present context. 

A phenomenological symmetry analysis of possible 
orderings in Ri-^Ca^MnOa was carried out in Rcf. 
[r?| . The phases that we find fit exactly into the phe- 
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nomenological classification that these authors propose. 
They conclude that three phases may exist: one with 
P2!/m symmetry, which coincides with our SCDW CE- 
structure; another of the Pmn2! type, corresponding 
to our BCDW orthogonal (Zener polaron) state, and a 
structure with Pm-symmetry that corresponds to our fer- 
roelectric intermediate phase. Experimental structural 
refinements do show that the actual symmetry of the low- 
temperature phase of Pr .6Ca .4MnO3 is Pm 0, i.e. it is 
in the intermediate phase. Our calculations provide the 
microscopic basis for this phenomenological treatment 
and shows that, indeed, all three phases are stable in a 
doping region close to x ~ 0.4. Apparently, with increas- 
ing x Pri-^CazMnOa goes from the _L or intermediate 
phase to the CE-phase at around x = 1/2, and then to 
the C-phase for overdoped samples. The (LaCa)Mn03- 
system is, due to a larger hopping t, apparently below the 
_L phase in Fig. 2: with increasing doping it goes from F 
to CE and then to the C-phase. In our phase diagram 
the A-phase appears between the metallic F and " insulat- 
ing" _L and CE-phases 0, whereas experimentally the 
A-phase usually exists only for x > 0.5 0|; the reason 
for this discrepancy is not clear at present. 

We should remark that we only considered different 
homogeneous states. It is known, however, that in many 
correlated electron models for the manganites the doped 
systems may be unstable with respect to phase separa- 
tion (PS) |20(. Also in our analysis we find for the in- 
termediate phase a convex dependence of the energy on 
the electron density, i.e. there is the tendency for PS 
into a pure CE and a pure _L phase with their appro- 
priate electron densities. However, the inclusion of the 
long-range Coulomb interaction in model calculations of 
the electronic properties prevents large-scale PS. Such a 
long-range interaction has to be taken into account in 
the case of PS, in order to guarantee the total charge 
neutrality of the system. These Coulomb interactions al- 
low only for the formation of weak inhomogeneities with 
smooth boundaries, which implies that the properties of 
the resulting PS state are similar to the properties of the 
homogeneous state. Nevertheless, the question of what 
the detailed properties of such an inhomogeneous state 
would be, is an interesting problem and deserves further 
investigation. 

Conclusions. We have shown that in manganites with 
doping concentrations close to i = 1/2 a novel type 
of magnetic order can appear besides the well known 
CE-type ordering. In this (_L) phase the neighboring 
spins are, depending on the direction of the bond, ei- 
ther oriented parallel, antiparallel or perpendicular to 
each other. Due to double exchange and orbital order- 
ing, parallel spins form strong dimer-like bonds, which 
we identify as the Zener polaron state that was recently 
observed experimentally close to x = 1/2 |4J. In consid- 



erable part of the phase diagram an intermediate phase, 
which is a superposition of a CE-type and orthogonal 
state, is stable. We predict that in this charge ordered 
intermediate state inversion symmetry is spontaneously 
broken and ferroelectric moments appear. This may be 
the first example of the appearance of ferroelectricity due 
to charge ordering in transition metal oxides. 
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